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ABSTRACT 

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of 
the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of 
Herschel, we examine the dust temperature and dust mass surface density distribution. We 
find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses 
compared to interarm regions. However, the distribution of dust temperature and mass are not 
spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the 
spiral arms lie ahead of the dust surface density peaks. The dust mass surface density corre- 
lates well with the distribution of molecular gas as traced by CO (J=3— >2) images (JCMT) 
and the star formation rate as traced by Ha with a correction for obscured star formation us- 
ing 24 /im emission. Using Hi images from THINGS to trace the atomic gas component, we 
make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean 
gas-to-dust ratio of 84 ± 4 with higher values in the inner region assuming a constant CO-to- 
H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor 
that varies with metallicity. 



1 INTRODUCTION 



Dust emission represents approximately one third of the bolomet- 
ric luminosity in normal spiral galaxies and dust acts to shape the 
spectral energy distribution (SED) of galaxies by attenuating short 
wavelength emission and re-emitting at longer wavelengths. Dust 
can also be used to trace the total gas content in galaxies, since dust 
and gas often coexist. An understanding of dust and its distribution 
are essential to reveal the underlying appearance of galaxies and 
the processes within (Draine 2003). 

Until recently it has been difficult to map the distribution and 
amount of dust in galaxies because, for most galaxies, only shorter 
wavelength infrared data has been available (i.e. X < 160 /im). 
Without longer wavelengths, only the warmer dust emission is 
measured and subsequently the dust mass is poorly constrained 
(Draine et al. 2007 and Galametz et al. 201 1). While some ground- 
based telescopes (i.e. SCUBA on JCMT; Holland et al. 1999 and 
SHARCII on the CSO; Dowel et al. 2003) have provided measure- 
ments in the far-infrared (350, 450 and 850 ^m) allowing for more 



accurate dust SED fits (e.g. Dunne et al. 2000, Dunne & Eales 2000 
and Dale et al. 2007), these images had low sensitivity due to back- 
ground noise from the atmosphere. 

The Herschel Space Observatory (Pilbratt et al. 2010) has al- 
lowed us to probe the colder dust component in galaxies using mea- 
surements at wavelengths spanning 70 to 500 fim at high sensitiv- 
ity and resolution. In comparison to the Spitzer Space Telescope 
(Werner et al. 2004), which also operated at 70 and 160 ^m (MIPS; 
Rieke et al. 2004), the point spread functions (PSFs) are smaller, the 
sensitivity is higher and there are fewer latent image effects. 

In this work we use Herschel observations to make dust tem- 
perature and mass maps for M83 which is part of the Very Nearby 
Galaxy Survey (VNGS; PI: C. Wilson). Due to its proximity (4.5 
± 0.3 Mpc; Thim et al. 2003), which allows for high spatial res- 
olution, and its relatively face-on orientation (inclination angle of 
24°), it makes an excellent candidate to study the spatial distribu- 
tion of the dust mass and temperature. 

M83 is a starburst galaxy with a prominent bar, which is out- 
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lined by dust lanes. Dynamical studies have shown that gas is fun- 
neled along the bar producing high rates of star formation at the 
centre (Knapen et al. 2010). M83 is rich in molecular gas with 13 
per cent of the disk mass consisting of molecular gas (Lundgren 
et al. 2004). Studies of the spiral arms have shown that tracers of 
the gas and star formation rate (SFR) are often offset from each 
other suggesting that the spiral arms dynamically induce star for- 
mation (Lord & Kenney 1991). However, the situation is compli- 
cated in that molecular clouds are also seen coincident with young 
star forming regions (Rand et al. 1999). Studies of the mid-infrared 
emission, which traces the very small dust grains and PAHs, have 
revealed a tight correlation with the radio continuum, which may 
suggest the anchoring of magnetic fields in the photoionizing shells 
of clouds (Vogler et al. 2005). The morphology of M83 and the fact 
that it has coverage at many other wavelengths allows for compar- 
isons of the varying environments within the galaxy and compar- 
isons of the dust distribution to tracers of the gas and SFR. 

We estimate the dust temperature and mass of M83 on angu- 
lar scales of 36 arcsec (790 pc) using 70 fim and 160 fim images 
from the Photodetector Array Camera and Spectrometer (PACS; 
Poglitsch et al. 2010) and the 250, 350, 500 /im maps from the 
Spectral and Photometric Imaging REceiver (SPIRE; Griffin et al. 
2010). We compare the spatial distribution of the dust temperature 
and mass with tracers of the atomic gas, molecular gas and star 
formation rate for different regions in the galaxy including the nu- 
cleus, bar and spiral arms. In §2 we describe how the infrared, gas 
and SFR images were rendered and processed. Our derivation of 
the dust temperature and dust mass distributions through SED fit- 
ting on each pixel is described in §3 and compared to the distribu- 
tions of molecular and total gas as well as the star formation rate 
surface density in §4. Our conclusions are summarized in §5. 



2 OBSERVATIONS & IMAGE RENDERING 
2.1 Herschel Images 

We use PACS images that are processed using both HIPE v5 and 
SCANAMORPHOS vif] (Roussel et al., sub.) and SPIRE images pro- 
cessed with HIPE and BriGAde (Smith et al., in prep). Details on 
how the images are processed can be found in Bendo et al. (201 1), 
the SPIRE Observer's Manua0and the PACS Observer's ManuaQ 
We briefly describe the steps here. The PACS images were cor- 
rected from the v5 photometric calibration files to v6 (with correc- 
tive factors of 1.119 and 1.174 for the 70 and 160 maps respec- 
tively) and were converted to Jy/sr from Jy/pixel . The SPIRE im- 
ages are first converted to Jy/pixel from Jy/beam using the respec- 
tive beam sizes (423 ± 3, 751 ± 4, 1587 ± 9 arcsec 2 , for the 250, 
350 and 500 /im images respectively). The images were then con- 
verted to Jy/sr. The SPIRE images were multiplied by the following 
values to convert from monochromatic intensities of point sources 
to monochromatic extended sources, 0.9939/1.0113, 0.9898/1.0060 
and 0.9773/1.0065 for the 250, 350 and 500 /im maps respectively. 
Table 1 lists the pixel sizes and FWHM of the PSFs for both the 
PACS and SPIRE images. Sky subtractions for all five images were 
made by subtracting the median value of several off-galaxy aper- 
tures. Fig.[Tjshows the five images in their native resolution, before 
convolution and deprojection. 



Table 1. Characteristics of Images 



Image 


Pixel Size [ arcsec] 


FWHM of PSF[ arcsec] 


70 /im 


1.4 


6.0 


160 fim 


2.85 


12.0 


250 fim 


6 


18.2 


350 fim 


8 


24.5 


500 fim 


12 


36.0 



Table 2. Adopted Properties of M83 



Inclination 


24° 


Position Angle 


225° 


Distance 


4.5 Mpc 



The images are all aligned and convolved to the 500 fim im- 
ages using kernels developed by Aniano et al. (201 1) for PACS and 
SPIRE. The final pixel size is 12 arcsec. However, in our quan- 
titative analysis we use 36 arcsec pixels, which correspond to the 
FWHM of the 500 fim map. We also deproject the images to a face- 
on orientation adopting the inclination and position angles from 
Tilanus & Allen (1993) (see Table 2). 

The uncertainties for each image are a combination of the un- 
certainty image files, which are corrected for convolution and de- 
projection, sky background uncertainties and the calibration uncer- 
tainties (0.03, 0.05, 0.07, 0.07, 0.07 for the calibration uncertainties 
of the 70, 160, 250, 350 and 500 fim images respectively). 

While Multiband Imaging Photometer images from Spitzer 
(MIPS; Rieke et al. 2004) are available for M83 at 70 fim and 
160 /im, they suffer from latent image effects and saturation in 
the nucleus due to M83's brightness (see Bendo et al. 2011 for a 
discussion). Thus, we choose not to include these images in our 
analysis. 

2.2 HI, CO & Gas Maps 

We use 21 cm line emission maps from THINGS (The HI Nearby 
Galaxy Survey; Walter et al. 2008) to trace the atomic gas. We use 
the natural weighted maps with a resolution of 13 arcsec and con- 
vert the integrated intensity to a surface density following Leroy 
et al. (2008). The HI maps are aligned and convolved to match the 
500 /im images and are deprojected according to the values in Table 
2. 

Muraoka et al. (2007) mapped CO (J=1^0) out to 3.5 
kpc. However, archival CO (J=3— s>2) data from the James Clerk 
Maxwell Telescope (JCMT; project: M06BN05) with a resolution 
of 14.5 arcsec allow us to map molecular gas distribution out to 5 
kpc. Thus, we use the CO (J=3— >2) but first scale the images to 
match the CO (J=l — s>0) of M83 using the ratio found by Muraoka 
et al. (2007), which varies with radius. Muraoka et al. (2007) found 
that the ratio of CO (J=3->2) to CO (J=l->0) intensity is as high as 
1.0 in the very inner regions and decreases to 0.65 at 0.5 kpc. Be- 
yond 0.5 kpc the ratio is roughly constant at 0.65. We use a radially 
varying ratio, S. At each radial position, R, between and 0.5 kpc 
we have 



1 SCANAMORPHOS Documentation: http: //www2 . iap. fr/users/roussel/herschel/ S=l.0 — 0.7R R < 0.5. (1) 

2 SPIRE Observer's Manual: [http : / /herschel . esac. esa. int/Docs/SPIRE/html/spire_om. html| 

3 PACS Observer's Manual: http : //herschel . esac . esa . int /Docs /P^Bey/OndnQ.^pkpC3 WfinusetSi0.65. Details on how the CO (J=3— >2) 
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Figure 1. M83 at 70 fim and 160 fim imaged using PACS (upper panels from left to right) and at 250 fim, 350 fim and 500 (im imaged using SPIRE (lower 
panel from left to right). The images are in Jy/sr and are scaled logarithmically with contours at equal increments (6, 6.5, 7, 7.5, 8, 8.5 and 9). The size of the 
PSF is illustrated in the top right hand corner of each image (see Table 1 for pixel and PSF sizes). 



images were processed can be found in Warren et al. (2010). We 
convert from antenna temperature to main beam temperature using 
a factor t\mb = 0.6. To estimate the surface density of H2 from 
CO we then use a constant conversion factor of X=2.0x 10 mol 
cm' 2 (K km s -1 ) -1 which corresponds to 3.2 Mq pc~ 2 (K km 
s ) (Wilson et al. 2009). We examine a conversion factor that 
varies with metallicity in §3.3. The H2 maps are aligned and con- 
volved to match the 500 fim images and are deprojected according 
to the values in Table 2. 

Together the HI and H2 maps are combined to make total gas 
surface density maps using a factor of 1.36 to account for helium. 



2.3 SFR Maps 

Measuring the SFR requires accounting for both unobscured star 
formation and dust obscured star formation. In this work we use 
Ha images, which are corrected for dust obscuration in compact 
sources using 24 fim emission maps. We use continuum subtracted 
Ha images from SINGG (Survey for Ionization in Neutral Gas 
Galaxies; Meurer et al. 2006) and 24 fim observations from En- 
gelbracht et al. (2005) that are reprocessed by Bendo et al. (2011, 
in prep.). The 24 fim emission traces star formation in compact re- 
gions that are obscured by dust, but also has a diffuse component 
from post-main sequence stars. Here we exclude the diffuse compo- 
nent and only including the 24fim emission from compact sources. 
Our goal is to trace star forming regions rather than radiation as- 



sociated with star formation. In Appendix B we compare our star 
formation tracer to linear combinations of images, which include 
the diffuse component. 

Several authors have demonstrated that 24 fim emission 
from unresolved s ources is correlated with the Ha emission from 
these sources (e.g.lCalzetti et al.ll2005ll2007l;|Pre"scott et alj[2007h . 
ICalzetti et alj d2007 ) derived the following equation to use 24 /im 
flux density measurements to correct the Ha flux from HII regions 
within a galaxy to produce an extinction-corrected Ha flux: 



/if a, corrected — /ffa observed + (0.031 ± 0.006) /24pn 



(2) 



We use the extinction-corrected Ha fluxes given by this equation 
as a star formation tracer in our analysis. However, to create a map 
of the extinction-corrected Ha emission, we first identify the un- 
resolved Ha sources and their 24 fim counterparts. We then apply 
the correction to each source and re-map the extinction-corrected 
emission into a new map. 

To identify the unresolved Ha so urces, we used SExtractor 
(version 2.5.0; iBertin & Arnoutslll996l) in double-image mode to 
create the initial source catalogues where detection and setting of 
apertures is done in the Ha image. All Ha images were first con- 
volved to the resolution of the 24 ^m image using a customised 
kernel created using the procedure outlined in Bendo et al. (201 1). 
The images were then regridded to the pixel size of the 24 fim map. 
Sources were selected when at least 20 pixels were 3a above the 
noise level in the Ha image used for detection. All photometry was 
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measured in apertures with a 40 arcsec (~27 pixels) diameter. We 
used this fixed radius because it is larger than the peak in the second 
Airy ring of the 24 data and the aperture corrections should be 
less than or equal to the calibration uncertainties. To minimize the 
number of faint spurious sources in our catalog, we exclude regions 
beyond the galaxy optical disk and foreground stars that appear as 
sources with positive and negative counterparts in the Ha images. 
We find that when run SExtractor is run in single-image mode, the 
mean absolutde value difference between the 24 /im flux densities 
using the two methods is ~13%. Thus, we add a 13% systematic 
error to all of the 24 measurements. 

After performing these measurements, we corrected the Ha 
measurements for foreground dust extinction within the Milky Way 
using A R = 0.176, which was calculate d by the NASA/IPAC 
Extragalactic Database based on data from Schl egel et ail (l998). 
We then applied Equation [2] to calculate extinction-corrected Ha 
fluxes, remapped the Ha fluxes into individual pixels in a map with 
very small (1 arcsec) pixels, and then convolved the data with a 
Gaussian function with a FWHM of 36 arcsec to match the PSF of 
the 500 /im map for analysis. 

Using the calibration of Kennicutt (1998) and Calzetti et al. 
(2007), we convert the corrected Ha maps, which we label as 
Ha jCor into an SFR: 

SFR(M yr" 1 ) = 5.3 x 10" 42 L corr (Ha, cor ) (erg s" 1 ). (3) 

For the remainder of this text we will refer to the SFR using this 
description. In Appendix B we also compare this tracer to other 
tracers derived from a linear combination of images. We examine 
an SFR made from the combination of the far-ultraviolet (UV) and 
24 fim emission as well as a combination of Ha and 24 fim emis- 
sion across the whole image without the compact source identifica- 
tion outlined above. 

All of the gas, dust and SFR maps are rebinned into 36 arcsec 
pixels corresponding to a physical size of 790 pc for the purposes 
of quantitative analysis (see Bendo et al, 201 1 for a discussion). 36 
arcsec represents the FWHM of the convolved data (see Table 1). 



3 DUST MAPS 

3.1 Temperature and Emissivity 

The dust temperature is determined by fitting each pixel in the 
PACS and SPIRE images with a single modified blackbody func- 
tion of the form: 

S\ = NBx(T)\- p , (4) 

where S\ is the flux density, B\(T) is the Planck function, N is a 
constant related to the column density of the material and f3 is the 
emissivity. We fit the points using MPFIT, a least-squares curve fit- 
ting routine for IDL (Markwardt 2009). Because the dust emission 
does not have a flat spectrum across the bands, during our fitting 
routine, we apply colour corrections to the five images. The appro- 
priate colour corrections were determined by first fitting the SED 
of the five images and determining the spectral index or slope at 
each point. The photometric colour corrections were determined 
by interpolation from the tables provided in the observer manuals 
of PACS and SPIRE. This process was repeated iteratively. 

The grain emissivity parameter (/}) of the dust reflects its prop- 
erties and is important for determining the dust temperature and 
mass. Typically, it has a value of ~2 (Draine & Lee 1984). How- 
ever, it can vary from as low as 1 to as high as 3 (Shetty et al. 2009a 



and references therein). While constant values of f3 of either 1.5 or 
2 have been used in many studies, /3 and the dust temperature may 
be anti-correlated due to the properties of amorphous grains (Dupac 
et al. 2001; Meny et al. 2007; Paradis et al. 2010). Thus, fixing /3 
can lead to erroneous temperature values. However, noise and line- 
of-site temperature variations can also produce this inverse relation 
making uncovering the intrinsic relation between j3 and tempera- 
ture challenging (Shetty et al. 2009a; Shetty et al. 2009b; Malinen 
et al. 2010). 

In this work we fit the dust emission using both a constant and 
a variable value of /3. For a variable /3 we permit values between 
1.0 (the limit of the Kramers-Kronig dispersion relation) and 2.5. 
The mean value of /3 over the entire galaxy is 2.1 ±0.1. For our 
constant /? fit, we assume a value of 2.0, which is close to the mean 
value found when /3 was allowed to vary. The mean \ 2 over all the 
fit pixels when using a variable f3 is 0.5 with a standard deviation 
of 0.3. In the case of a constant j3 the mean value is 0.7 with a 
standard deviation of 0.3. These \ 2 values are small suggesting 
that our errors may be large. This is likely due to the calibration 
errors, which are correlated. Fig.[2]shows four sample fits in log-log 
space for pixels in the nuclear, bar, spiral arm and interarm regions. 
The curves are all normalized to the peak value. The bottom panel 
shows how the data deviates from the fit (i.e. log S\ — log 7a). In 
each panel we list the fitted temperature and emissivity. 

Fig. |3] shows the spatial distribution of f3 when it was allowed 
to vary (left) and the resulting temperature maps for a variable /3 
(centre) and constant, /3 = 2 (right). The high /3 values on the edge 
of the map are likely due to the low signal-to-noise in these regions. 
One notes clear differences in the two temperature maps, particu- 
larly in the inner regions. With a constant f3, the central temperature 
peak is elongated perpendicular to the bar and shows an outer ring- 
like structure, whereas with a variable /?, the central peak is much 
more symmetrical. Over much of the galaxy, the value of f3 was 
close to 2, but in the inner regions there is a ring of higher f3 values. 
This ring is the source of the differences in the inner regions of the 
two temperature maps. Bendo et al. (2003) also found values of /3 
close to 2 in the inner regions (within 135 arcsec) of M83 using ISO 
(Kessler et al. 1996) and SCUBA (Holland et al. 1999) data. M83 
had the steepest emissivity law in their analysis of eight galaxies 
using a variable j3. 

Fig.[4]shows h° w /? varies with temperature. There is an anti- 
correlation between /? and temperature. This is in agreement with 
other studies (i.e. Dupac et al. 2001, Paradis et al. 2010 and Dariush 
et al. 201 1) and recent results from Planck, which have revealed an 
anti-correlation between /3 and temperature from SED fits of the 
integrated emission of nearby galaxies (Planck Collaboration et al. 
2011). 

However, whether we have uncovered the intrinsic inverse re- 
lation remains unclear as noise and line-of-sight temperature vari- 
ations can also produce a similar relation (Shetty et al. 2009a and 
2009b). While there is a trend of decreasing f3 with increasing tem- 
perature, we note that in the inner regions (dark blue points) where 
there are mostly high temperatures (T > 25 K), there are some re- 
gions with j3 exceeding 2. These values lie in the ring-like structure 
seen in the left panel of Fig.[3]and deviate from the anti-correlation. 
While it is possible that the ring is simply an artifact of imperfect 
PSF matching, it may also represent a region of different physical 
properties. 

Using a single temperature modified blackbody fit means that 
at each pixel we have an average dust temperature. Typically, galax- 
ies have a mixture of dust at different temperatures with a warmer 
and cooler component. Here we briefly examine to what extent 
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these features might affect our results. In order to examine the 
possibility of a warmer component, we fit the SED with our con- 
stant /3 value of 2 dropping the 70 fim point. If the 70 fim point is 
higher than our fit using 160, 250, 350 and 500^im, this suggests a 
warmer component (i.e. Bendo et al, 2010a and Smith et al. 2010). 
We find that the SED fit over-predicts the 70 value over much 
of the galaxy. This suggests that including the 70 fim point is nec- 
essary in order to properly constrain the peak. By examining the 
far-infrared colour variations in M83, Bendo et al. (2011) found 
that the dust could be divided into a component radiating primarily 
at 250-500 and another at 70-160 fim, which is not evident in 
the test we perform here. 

There may also be a cold dust component (T~5-10K), which 
has been attributed to a so-called "submm excess" (i.e. Galliano et 
al. 2005, Galliano et al. 2005, Bendo et al. 2006, Zhu et al. 2009, 
O'Halloran et al. 2010 and Gordon et al. 2010). This submm ex- 
cess has also been explained, in the case of the Milky Way, by 
wavelength-dependent variations in the dust emissivity (i.e. Reach 
et al. 1995 and Paradis et al. 2009). Typically, one requires mea- 
surements at longer wavelengths (i.e. 850 mum) in order to iden- 
tify the excess, but even at 500 fim the excess may present itself 
(i.e. Gordon et al. 2010). We repeat a similar procedure to that out- 
lined above for the 70 fim observations. In general, we find that the 
500 fim observed value exceeds the predicted value of the SED fit 
using the other four wavelengths by ~10 per cent. This is in agree- 
ment with the excess found by Gordon et al. (2010). Fig.[5]shows 
the percent difference between 500 fim (right) point and the fit us- 
ing the other four maps. Positive values indicate an excess. We note 
that the regions where the 500 fim is under-predicted correspond to 
the ring-like structure in the variable f3 map. Since we were forced 
to keep /3 fix in order to properly fit the SED with only four points, 
this ring may be an artifact of fixing j3. We stress that the excess 
found here is small in comparison to our uncertainties and longer 
wavelength measurements are necessary to confirm this finding. 

3.2 Dust Mass 

We calculate the dust mass using the following equation at 250 fim: 

S250D 2 



Mdust — 



(5) 



K250B(\250,T) 

where S\ is the flux density from our SED fit at our chosen 250 fim 
wavelength (i.e. not the value in the 250 fim map), D is the dis- 
tance to M83 (4.5 Mpc; Thim et al. 2003), B(\, T) is the Planck 
function using 250 and the best fit temperature and k is the 
assumed emissivity. We use a k of 2.72xl0 5 (\/fim)~ 2 cm 2 g -1 
where A=250 fim. The value of n is taken from Li & Draine (2001) 
and modified with a correction discussed in Draine (2003). 

We calculate the mass at 250 as it is not overly sensitive to 
the temperature like the shorter PACS wavelengths and has lower 
uncertainties than the other SPIRE wavelengths. We tested the other 
wavelengths as fiducial wavelengths for the dust mass maps and 
found that the total dust masses agree within the uncertainties. 

We compared the spatial distribution of the dust surface den- 
sity using either a constant /3 or variable f3 and found that they 
are very similar. We find a total dust mass of 4.0 ± 0.6 x 10 7 Mq 
with a variable j3 and 3.9 ± 0.6 x 10 7 M Q with a constant j3 within 
r=300 arcsec. The difference between the two dust mass maps in in- 
dividual pixels is at the 10 per cent level. Galametz et al. (201 1) did 
global SED fits of M83 using infrared data from IRAS and submm 
data at 540 fim from Hildebrand et al. 1977. Without the submm 
data, they found a total dust mass of 8.5 x 10 8 Mq . However, when 
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Figure 2. Sample modified blackbody fits for regions in the nucleus, bar, 
arm and interarm with the colour-corrected PACS and SPIRE points using a 
variable (3 (see centre of each panel for fit value). The peak of the function 
determines the temperature (listed in each panel). The functions are shown 
in log-log and are normalized to the peak of the function. The uncertainties 
are smaller than the symbol sizes. The deviation between each point and the 
fit (i.e. log S\ — log I\) is shown below each panel. 



they included the submm data, the total dust mass decreased to 8.5 
x 10 6 Mq, illustrating the need for longer wavelength data in or- 
der to properly constrain the SED fit. Our dust mass, while higher, 
is much closer to their global fit using the older longer wavelength 
data. We also compared the total dust mass from the fit, without 
using the 70 fim and 500 fim maps. In these cases, we assumed 
a j3 value of 2. In both cases, the total dust mass was identical to 
that using all five maps within the uncertainties. In Appendix A we 
also compare our dust mass estimations with those found using the 
models of Draine & Li (2007). 



4 RESULTS & DISCUSSION 

4.1 Distribution of Dust Temperature and Mass 

Figs.[3]and Fig.[6]show the dust temperature and dust mass surface 
density maps respectively for M83 based on the modified black- 
body fits with a variable /3. The highest temperatures (~ 30K) are 
found in the central peak, the spiral arms and the southwestern bar- 
spiral arm transition. This is likely due to the increased star forma- 
tion activity in these regions and, in the case of the bar end, orbit 
crowding of the molecular gas (Kenney & Lord 1991) and subse- 
quent star formation activity. Regions outside of the nucleus, but 
along the bar show lower temperatures, which are comparable to 
the temperatures in the interarm regions. 

The dust mass surface density is shown in Fig. [6] in log solar 
masses per square parsec. The regions of highest dust concentra- 
tion are located in the nucleus and at the ends of the bar. The right 
panel shows the dust surface density with the temperature contours 
overlaid. The peak temperature in the south-western bar-spiral arm 
transition peak and the north-eastern spiral arm are offset from the 
dust mass surface density peaks. The dust temperature is highest 
ahead of the dust mass surface density arms if we assume a trailing 
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Figure 3. Spatial distribution of /3 (left), dust temperature using a vai'iable ft in the modified blackbody fit (centre) and dust temperature using a constant of 
2 in the modified blackbody fit (right). All three images show the same field of view. For scale we show a 1' distance on the deprojected images (left panel). 
The contours in the /3 map are at 1.5, 1.8, 2.1 and 2.4. The contours in the two temperature maps (centre and right panels) are the same and lie between 22 and 
30 K in increments of 2. Over much of the galaxy the variable ft is roughly 2. However, there is an inner ring of higher ft values which leads to the differences 
in the central regions of the two temperature maps. 
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Figure 4. ft versus temperature for the variable ft SED fits. Each point rep- 
resents a 36 arcsec region corresponding to the FWHM of the PSF in the 
convolved maps. The colours mark the radial position of the points with 
blue and green points within 3 kpc of the centre and yellow and red points 
extending further out. The black point illustrates the mean uncertainties of ft 
and the temperature. The solid line shows the two-degree polynomial curve 
fit to the points, taking into consideration the uncertainties. We find that ft 
and the temperature are anti-correlated in agreement with other studies. It is 
unclear whether this anti-correlation is intrinsic or is simply a reflection of 
noise and line-of-sight temperature variations (Shetty et al. 2009a and b). 



spiral structure rotating in the clockwise direction. Dust tempera- 
ture and mass maps for NGC 4501 and NGC 4567/8 in the Herschel 
Virgo Cluster Survey have also shown that the peaks in dust tem- 
perature were not coincident with those in mass (Smith et al. 2010). 



Figure 5. The percent difference between the flux of the 500 (im map and 
that predicted by the fit using the 70, 160, 250 and 350 fim maps. Positive 
values mean the observed value is greater than predicted, while negative 
values mean the observed is less than predicted. With the exception of the 
inner ring, much of the galaxy shows an excess of 500 /im emission. 



4.2 Comparison with Gas & SFR 

We now compare the spatial distribution of the dust mass surface 
density and temperature with tracers of the atomic and molecular 
gas and the SFR. In Fig.[7]we show the tracers with contours of the 
dust surface density (upper panel) and contours of the dust temper- 
ature (lower panel). In Fig. [8] we plot the pixels in the dust mass 
surface density (upper row) and temperature maps (lower row) ver- 
sus those in the different tracers. Each point represents an area of 
36 arcsec (the FWHM of convolved beam) and we denote regions 
in the nucleus, bar and spiral arms in orange, purple and cyan re- 
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Figure 6. Dust mass surface density map in log Mq pc -2 for M83 with contours (left) and with contours of the temperature map overlaid (right). The 
temperature map is that made with a variable /3 (see middle panel of Fig. [5}. In many cases, peaks in temperature are offset ahead of peaks in mass in the spiral 
and bar structures. 



spectively. The choice for these regions was made by-eye using 
3.6 fim images from IRAC (Spitzer Local Volume Legacy; Dale 
et al. 2009) that were convolved to the 500 fim resolution and de- 
projected. 3.6 maps roughly trace the underlying stellar mass 
density, which makes them well-suited for probing morphological 
structures. Although PAH and hot dust emission can also contribute 
in this band in regions dominated by emission associated with star 
formation (Mentuch et al. 2010), since these regions in M83 al- 
ready correspond to regions with enhanced stellar mass densities, 
this contamination does not affect our qualitative approach in de- 
composing the galaxy's morphology. 

Examining the upper panel of Fig.|7jwe see that peaks in the 
dust mass surface density are virtually coincident with those of the 
molecular gas and the SFR. The position of the spiral arm structure 
in the dust mass surface density is also coincident with that in the 
atomic gas map. In the nucleus, however, we find little atomic gas. 
In contrast to the excellent spatial correspondence between the gas 
and SFR tracers with the dust mass surface density, we note that 
in the lower panel of Fig. [7] the dust temperature does not have as 
close of an agreement. 

In Fig. [8] we examine in a more quantitative fashion the corre- 
lations between the dust and the different gas and SFR tracers (see 
the correlation coefficients in the upper right of each panel). Again 
here we find that the molecular gas is tightly correlated with the 
dust mass surface density. While the total gas is also tightly corre- 
lated with the dust surface density, it is important to note, that over 
the radii considered, the molecular component dominates the mass 
budget (Crosthwaite et al. 2002) so this is not too surprising. We 
note that the dust temperature tends to increase with rising molecu- 
lar gas and SFR, but remains flat with the atomic gas density. This 
is simply the result of the increasing atomic gas density with ra- 
dius, but decreasing dust temperature with radius (see Fig. [9] for 
radial profiles). The dust temperature is somewhat correlated with 
the molecular gas and SFR, but the correlations are not as strong 
as with the dust mass surface density. Bendo et al. (201 1) has also 
found offsets between colour temperatures (160/250 and 250/350) 



with the Ha spiral arm structure. This may reveal that while the 
dust mass surface density is tightly correlated to molecular clouds 
and thus, star formation, dust heating may be caused by a variety 
of sources and not just recent star formation activity. Stars form in 
dense, cold molecular clouds. After the stars have formed, they heat 
the surrounding dust that lies within one mean free path. Thus, it is 
possible that colder dust lies in the star forming regions. As the dust 
disperses it will be heated by more stars in the surrounding region 
and will subsequently become hotter. 

Examining the different regions of the galaxy reflected by the 
colour-coding of the points, we find that the nuclear regions tend to 
have higher dust temperatures and higher dust masses for the same 
value of the gas surface density. The latter, however, may be due to 
radial variations in the gas-to-dust ratio related to metallicity (see 
§4.5 for a further discussion of this). 



4.3 Radial Profiles 

Fig.|9]shows the radial profiles of the gas tracers, SFR, dust temper- 
ature, dust mass surface density and dust emissivity profiles. Here 
again we mark the pixels in the nuclear (orange), bar (purple) and 
arm (cyan) regions based on the 3.6 /im map. The solid line de- 
notes the azimuthal average using radial annuli of 12 arcsec width. 
The dashed lines show the standard deviation around the azimuthal 
average. The points represent 36 arcsec pixels. 

The molecular gas profile of M83 (centre panel of Fig. [9]l is 
typical of a strongly barred spiral galaxy. The profile shows a cen- 
tral peak of molecular gas that declines, but then rises again pro- 
ducing a "shoulder". The shoulder forms at the bar-spiral arm tran- 
sition region (Nishiyama et al. 2001) and is also seen in the SFR, 
dust temperature and dust mass surface density profiles. In inves- 
tigating the UV extinction profiles of M83, Boissier et al. (2005) 
also found that the SFR and extinction rates were enhanced at the 
bar-spiral arm transition region. 

We identified the peak of the shoulder in the different profiles 
by finding the first maximum beyond 1 kpc. It is marked in each 
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Figure 7. Atomic (left), molecular (second from left) and total (third from left) gas surface density (logarithmic scale) in Mq pc — 2 and the SFR (right) 
as determined from dust-corrected Ho emission (as described in §2.3) in Mq kpc~ 2 yr _1 . Contours of the dust mass surface density (upper panel) and 
temperature (bottom panel) distribution are overlaid. We find close agreement between the dust mass surface density and molecular gas and SFR distribution. 
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Figure 8. The dust temperature (lower panel) and dust mass surface density (upper panel) versus atomic gas (left), molecular gas (second from left), total gas 
(third from left) mass and the SFR (right) as traced by dust-corrected Ha emission. Each point represents a pixel 36 arcsec in size (the FWHM of the convolved 
maps). The pixels in the nucleus, bar and spiral arms are denoted in orange, purple and cyan respectively. The black pixels are in the interarm region or are 
beyond the extent of the clearly defined spiral structure. The correlation coefficients between the different tracers are listed in the upper right of each panel. The 
dust mass surface density is tightly correlated with the molecular gas, total gas and SFR. The dust temperature is surprisingly only somewhat correlated with 
the SFR. While a rise in temperature is associated with increased molecular gas and SFR, the atomic gas mass is associated with a decrease in temperature. 
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Table 3. Shoulder Position 



Profile 


Shoulder Position [kpc] 


Dust Temperature 


2.9 


Dust Mass 


2.1 


SFR 


2.4 


Molecular Gas 


2.1 



plot by the dashed vertical line and it can also be seen in the tran- 
sition from purple to cyan squares. The shoulder does not seem to 
be associated with the position of corotation, which is beyond the 
bar end for M83 and is marked in these plots with the solid vertical 
line at 3.7 kpc (Lundgren et al. 2004). 

The shoulder is associated with high dust temperature and 
mass surface densities in the bar-spiral arm transition (especially 
the southwestern arm) as seen in the spatial distribution maps in 
§4.1. The position of the shoulder varies with the different profiles 
(see Table 3 for the positions). The shoulder in the molecular gas 
and dust mass surface density profile are coincident. As we dis- 
cussed in §4.3, there is close spatial agreement between the two 
maps. Both the dust temperature and SFR have shoulders that are 
further out in the disk, and surprisingly they are not coincident. 

We note in the bottom right panel Fig. [9] that the tempera- 
ture of the dust in the bar region right before the shoulder is much 
higher than average (denoted by solid curve). This is likely due to 
increased SF due to orbit crowding in this region. 

The top right panel of Fig. [9] shows how the value of ft varies 
with radius. /3 is close to 2 over much of the disk, with higher values 
in the inner 1 kpc. 



4.4 Gas-to-Dust Ratio 

Fig. [10] shows the radial variation of the gas-to-dust ratio (GDR), 
which is derived using the total gas mass (atomic and molecular) 
divided by the dust mass. Both the azimuthal average (black solid 
line) and the region within the standard deviation (grey shaded 
area) are shown. Across M83 we find a mean value of 84 ± 4 (stan- 
dard deviation: 40) which is roughly the galactic value of 100-200 
(e.g. Sodroski et al. 1997). The right vertical panel of Fig.fTOlshows 
the histogram of the GDR. Only a few regions in the galaxy have 
high GDRs in contradiction with previous estimates for M83 (De- 
vereux & Young 1990) which did not probe the cold dust. 

In the central regions, the GDR attains much higher values 
(over 200). This may simply be due to higher metallicities in the 
inner regions. One would expect a higher GDR with increasing 
metallicity. It may also be due to the Xco factor which can of- 
ten overestimate the molecular gas in the central regions of spiral 
galaxies including the Milky Way by a factor of 2 or more (i.e. 
Regan 2000). 

The Xco factor has been shown to be anti-correlated with 
metallicity in several studies (i.e. Wilson 1995, Arimoto et al. 1996, 
Boselli et al. 2002 and Israel 2005). A recent study by Magrini et al. 
(201 1) has shown that using an Xco factor that varies with metal- 
licity may have a large impact on the radial profile of the GDR. 
This is particularly true for the inner regions, where the molecular 
gas dominates. 

In light of this, we re-calculate our estimate of molecular hy- 
drogen using an Xco factor which varies with metallicity. Israel 
(2000) has derived a relationship between the Xco factor and the 
oxygen abundance as: 
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Figure 10. GDR versus radius (left) and histogram of the GDR (right). The 
azimuthal average of the GDR using a constant Xco factor (solid line) 
and a radially varying one (dashed line) are shown. The grey shaded region 
shows the area within one standard deviation of the azimuthal average. The 
histogram is normalized to the peak value and the mean GDR and standard 
deviation are listed (upper right). 



logXco = 12.2 



2.51ogg 



(6) 



We use the oxygen abundance gradients of Bresolin et al. 
(2009), which are based on linear regression fits of the results of 
several strong-line studies in the outer regions of M83 and com- 
piled studies of the inner region by Bresolin & Kennicutt (2002) 
and Bresolin et al. (2005). The abundance ratio decreases linearly 
from 12 + log 0/H=8.77 in the centre to 12 + log 0/H=8.6 at 6 
kpc. 

Fig. 1 101 shows the GDR using a variable Xco factor (dashed 
line). One notes that the GDR in the inner regions are much reduced 
and there is now a slight trend towards an increasing GDR with 
radius. Some of the values in the inner region are much lower (~60) 
than the galactic value of 100-200. This may not be too surprising 
given that the metallicity of M83 is super-solar in the inner regions 
(Bresolin et al. 2005), which also implies higher dust masses and 
subsequently a lower GDR. This shows the sensitivity of the GDR 
to an Xco factor that varies with metallicity. 

The consistency of the GDR over much of the disk beyond the 
inner regions is in contradiction to several other studies of similarly 
nearby spiral galaxies. Recent studies of M99 and Ml 00 (Pohlen et 
al. 2010) and NGC 2403 (Bendo et al. 2010b) have shown an in- 
creasing GDR with radius. Munoz-Mateo et al. (2009) also showed 
that almost all galaxies in the SINGS sample showed radially in- 
creasing GDRs. However, there were two exceptions. NGC 3627 
showed a decreasing GDR and NGC 2976 showed a fairly consis- 
tent GDR. 



5 CONCLUSIONS 

The close proximity of M83 has allowed us to probe in detail the 
dust temperature and mass distribution using images from PACS 
and SPIRE. With its prominent bar and spiral arms we have also 
been able to compare different environments and see the correspon- 
dence of the dust temperature and mass with tracers of the gas and 
SFR. We briefly summarize our findings: 

• We find an anti-correlation between /3 and temperature with 
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Figure 9. Radial profiles of the atomic gas (left, bottom), molecular gas (left, middle), SFR (left, top), dust temperature (right, bottom), dust mass surface 
density (right, middle) and emissivity, /3 values (right, top). The solid curve denotes the azimuthal average and the dotted lines show the standard deviation of 
this average. 36 arcsec pixels in the nuclear, bar and spiral regions are marked in orange, purple and cyan respectively. Corotation and the shoulder position 
(vertical solid and dashed line respectively) are marked in the molecular gas, SFR, dust temperature and mass surface density panels. 



a mean value of /3 of 2.1 ± 0.1. However, it is possible this anti- 
correlation is an artifact of noise and line-of-sight temperature vari- 
ations (Shetty et al. 2009a and 2009b). A comparison of the dust 
temperature maps made using a constant /3=2 and a variable f3 
shows that the most prominent differences are in the inner regions. 
The dust mass surface density maps are virtually identical and have 
the same total dust mass within the uncertainties; 

• We find a total dust mass of4.0±0.6 x 10 7 M Q ; 

• We find that the dust temperature and mass surface density dis- 
tributions are not spatially coincident. In particular, while the dust 
mass surface density correlates extremely well with the molecular 
gas and SFR, we find that the dust temperature peaks are offset. In 
the spiral arms, we find that the dust temperature peaks lie ahead of 
the dust mass surface density (assuming a trailing spiral structure); 

• The dust temperature is more strongly correlated with the 
molecular gas than the SFR; 

• The nuclear and bar-arm transition regions show enhanced 
dust mass surface densities and temperature. The bar-arm transi- 
tion regions lead to the development of a 'shoulder' in the radial 
profile plots. The position of the shoulder varies between the dust 
mass surface density and temperature radial profiles; 

• The GDR has a mean value of 84 ± 4, which is close to the 
galactic value. In the inner regions the GDR extends as high as 200; 

• Using an Xco factor that varies with metallicity decreases the 
GDR in the inner regions. 

K. Foyle acknowledges helpful conversations with Brent 
Groves, Frank Israel, Adam Leroy, Giovanni Natale, Karin Sand- 
strom and Ramin Skibba. This research was supported by grants 
from the Canadian Space Agency and the National Science and 
Engineering Research Council of Canada (PI: C. Wilson). PACS 
has been developed by a consortium of institutes led by MPE 



(Germany) and including UVIE (Austria); KU Leuven, CSL, 
IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAF- 
IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC (Spain). This de- 
velopment has been supported by the funding agencies BMVIT 
(Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR 
(Germany), ASI/INAF (Italy), and CICYT/MCYT (Spain). SPIRE 
has been developed by a consortium of institutes led by Cardiff 
University (UK) and including Univ. Lethbridge (Canada); NAOC 
(China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC 
(Spain); Stockholm Observatory (Sweden); STFC and UKSA 
(UK); Imperial College London, RAL, UCL-MSSL, UKATC, 
Univ. Sussex (UK); and Caltech, JPL, NHSC, Univ. Colorado 
(USA). This development has been supported by national fund- 
ing agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS 
(France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC 
(UK); and NASA (USA). HIPE is a joint development by the Her- 
schel Science Ground Segment Consortium, consisting of ESA, the 
NASA Herschel Science Center, and the HIFI, PACS and SPIRE 
consortia. This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the let Propulsion 
Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. 



REFERENCES 

Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. 2011, 
larXiv: 1106.50651 

Arimoto, N., Sofue, Y., & Tsujimoto, T. 1996, PASJ, 48, 275 
Bendo, G. L, et al. 2003, AJ, 125, 2361 

Bendo, G. J., Dale, D. A., Draine, B. T., et al. 2006, ApJ, 652, 283 
Bendo, G. J., et al. 2010a, A&A, 518, L65 



The Dust & Gas Properties ofM83 1 1 



Bendo, G. J., et al. 2010b, MNRAS, 402, 1409 
Bendo, G. J., Boselli, A., Dariush, A., et al. 2011, 
larXiv: 1109.02371 

Bertin E., Arnouts S., 1996, A&AS, 117, 393 
Boissier, S., Gil de Paz, A., Madore, B. R, et al. 2005, ApJL, 619, 
L83 

Boselli, A., Lequeux, J., & Gavazzi, G. 2002, A&A, 384, 33 
Bresolin, R, & Kennicutt, R. C, Jr. 2002, ApJ, 572, 838 
Bresolin, R, Schaerer, D., Gonzalez Delgado, R. M., & Stasinska, 

G. 2005, A&A, 441,981 
Bresolin, R, Ryan- Weber, E., Kennicutt, R. C, & Goddard, Q. 

2009, ApJ, 695, 580 
Calzetti D. et al., 2005, ApJ, 633, 871 
Calzetti D. et al., 2007, ApJ, 666, 870 

Crosthwaite, L. P., Turner, J. L., Buchholz, L., Ho, P. T. P., & 

Martin, R. N. 2002, AJ, 123, 1892 
Dale, D. A., Gil de Paz, A., Gordon, K. D., et al. 2007, ApJ, 655, 

863 

Dale, D. A., Cohen, S. A., Johnson, L. C, et al. 2009, ApJ, 703, 
517 

Dariush, A., et al. 2011 JarXiv:1106.6195l 
Devereux, N. A., & Young, J. S. 1990, ApJ, 359, 42 
Dowell, C. D., Allen, C. A., Babu, R. S., et al. 2003, Proc. SPIE, 
4855, 73 

Draine, B. T., et al. 2007, ApJ, 663, 866 
Draine, B. T., & Li, A. 2007, ApJ, 657, 810 
Draine, B. T. 2003, ARA&A, 41, 241 
Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 
Dupac, X„ et al. 2001, ApJ, 553, 604 

Dunne, L„ Eales, S., Edmunds, M., et al. 2000, MNRAS, 315, 115 
Dunne, L., & Eales, S. A. 2001, MNRAS, 327, 697 
Engelbracht, C. W., Gordon, K. D., Rieke, G. H., Werner, M. W., 

Dale, D. A., & Latter, W. B. 2005, ApJl, 628, L29 
Galametz, M., Madden, S. C, Galliano, R, et al. 201 1, A&A, 532, 

A56 

Galliano, R, Madden, S. C, Jones, A. P., et al. 2003, A&A, 407, 
159 

Galliano, R, Madden, S. C, Jones, A. P., Wilson, C. D., & 

Bernard, J.-P. 2005, A&A, 434, 867 
Gil de Paz, A., et al. 2007, ApJS, 173, 185 

Gordon, K. D., Galliano, R, Hony, S., et al. 2010, A&A, 518, L89 
Griffin, M. J., et al. 2010, A&A, 518, L3 
Holland, W. S., et al. 1999, MNRAS, 303, 659 
Israel, R 2000, Molecular Hydrogen in Space, 293 
Israel, F. P. 2005, A&A, 438, 855 
Kenney, J. D. P., & Lord, S. D. 1991, ApJ, 381, 118 
Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 
Kennicutt, R. C, Jr., et al. 2009, ApJ, 703, 1672 
Kessler, M. R, et al. 1996, A&A, 315, L27 
Knapen, J. H., Sharp, R. G., Ryder, S. D., et al. 2010, MNRAS, 
408, 797 

Leroy, A. K, Walter, R, Brinks, E., Bigiel, R, de Blok, W. J. G, 
Madore, B., & Thornley, M. D. 2008, AJ, 136, 2782 

Liu, G, Koda, J., Calzetti, D., Fukuhara, M., & Momose, R. 201 1, 
ApJ, 735, 63 

Lord, S. D., & Kenney, J. D. P. 1991, ApJ, 381, 130 

Lundgren, A. A., Olofsson, H., Wiklind, T., & Rydbeck, G. 2004, 

A&A, 422, 865 
Magrini, L., et al. 2011, larXiv:1106.0618l 

Malinen, J., Juvela, M., Collins, D. C, Lunttila, T., & Padoan, P. 
2011, A&A, 530, A101 



Markwardt, C. B. 2009, Astronomical Data Analysis Software 

and Systems XVIII, 41 1, 251 
Mentuch, E., Abraham, R. G, & Zibetti, S. 2010, ApJ, 725, 1971 
Meny, C, Gromov, V., Boudet, N., Bernard, J.-P, Paradis, D., & 

Nayral, C. 2007, A&A, 468, 171 
Meurer, G. R., et al. 2006, ApJS, 165, 307 
Munoz-Mateos, J. C, et al. 2009, ApJ, 701, 1965 
Muraoka, K, et al. 2007, PASJ, 59, 43 
Nishiyama, K., Nakai, N., & Kuno, N. 2001, PASJ, 53, 757 
O'Halloran, B., Galametz, M., Madden, S. C, et al. 2010, A&A, 

518, L58 

Paradis, D., Bernard, J.-P, & Meny, C. 2009, A&A, 506, 745 
Paradis, D., et al. 2010, A&A, 520, L8 
Pilbratt, G. L„ et al. 2010, A&A, 518, LI 

Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2011, 
larXiv: 1101.20451 

Poglitsch, A., et al. 2010, A&A, 518, L2 

Pohlen, M., et al. 2010, A&A, 518, L72 

Prescott M. K. M. et al., 2007, ApJ, 668, 182 

Rand, R. J., Lord, S. D., & Higdon, J. L. 1999, ApJ, 513, 720 

Reach, W. T., Dwek, E., Fixsen, D. J., et al. 1995, ApJ, 451, 188 

Regan, M. W. 2000, ApJ, 541, 142 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525 
Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 
154, 25 

Shetty, R., Kauffmann, J., Schnee, S., Goodman, A. A., & Er- 

colano, B. 2009a, ApJ, 696, 2234 
Shetty, R., Kauffmann, J., Schnee, S., & Goodman, A. A. 2009b, 

ApJ, 696, 676 
Smith, M. W. L., et al. 2010, A&A, 518, L51 
Sodroski, T. J., Odegard, N., Arendt, R. G., Dwek, E., Weiland, 

J. L., Hauser, M. G, & Kelsall, T. 1997, ApJ, 480, 173 
Thim, F, Tammann, G. A., Saha, A., Dolphin, A., Sandage, A., 

Tolstoy, E., & Labhardt, L. 2003, ApJ, 590, 256 
Tilanus, R. P. J., & Allen, R. J. 1993, A&A, 274, 707 
Vogler, A., Madden, S. C, Beck, R., et al. 2005, A&A, 441, 491 
Walter, F, Brinks, E., de Blok, W. J. G, Bigiel, F, Kennicutt, 

R. C, Jr., Thornley, M. D., & Leroy, A. 2008, AJ, 136, 2563 
Warren, B. E., et al. 2010, ApJ, 714, 571 
Wilson, C. D. 1995, ApJl, 448, L97 
Wilson, C. D., et al. 2009, ApJ, 693, 1736 
Werner, M. W., et al. 2004, ApJS, 154, 1 

Zhu, M., Papadopoulos, P. P., Xilouris, E. M., Kuno, N., & Lisen- 
feld, U. 2009, ApJ, 706, 941 



APPENDIX A: COMPARISON WITH DRAINE & LI 
MODELS 

We compare our dust mass estimation using the modified black- 
body analysis with that found using the models from Draine & Li 
(2007) and the five PACS and SPIRE images. The models are fit 
using the two-component emission model put forth by Draine & 
Li (2007). In their model, the mid- through far-infrared emission 
due to polycyclic aromatic hydrocarbons (PAHs), very small dust 
grains and carbonaceous dust grains is expressed as a function of 
the underlying interstellar radiation field (ISRF). The ISRF is rep- 
resented by two components, one consisting of an emission com- 
ponent modelling a photodissociation region (PDR) with a range of 
starlight intensities from U = U m i n to U=10 6 and an ISM compo- 



12 K. Foyle etal. 
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Figure Al. Dust mass surface density in log Mq pc -2 (left) and U m i n 
(right) based on fits using the Draine & Li (2007) models. While the dust 
mass surface density map looks very similar to that made with the modified 
blackbody fit, the surface densities are all approximately 1.5 times higher 
(see Fig.[6j. The U m ; n map shows similar features to the dust temperature 
map derived from a modified blackbody fit with a constant j3 (see Fig. [5}. 



of Ha and 24 fim (Calzetti et al. 2007). We first briefly describe the 
image rendering for this analysis. 

In order to make the SFR;yv r +24Mm, we convolved the UV and 
24 fim maps to the 500 fim resolution and deprojected them in 
accordance with Table 2. The UV and 24 images are combined 
in the following fashion: 

Xsfr = (8.1 x 10" 2 Ifuv + 3.2 x lO -3 ^), (Bl) 

where T,sfr has units of Mq kpc~ 2 yr _1 and the FUV and 24 fim 
intensity are each in MJy sr" 1 (for more details see Appendix A of 
Leroy et al. 2008). 

For the SFR£f a +24^m, we convolve each map to the 500 
and deproject them in accordance with Table 2. The Ha map is 
corrected for foreground dust extinction in the Milky Way as previ- 
ously described in §2.3. Rather than applying the 24 fim correction 
to HII regions within the galaxy, we simply combine the two maps 
using: 



nent with a single radiation field of U=Umir0. The bulk of dust 
mass is contained in this latter component. In our least-squares 
fitting procedure, Umin ranged from 0.1 to 25 and we vary the 
fraction, 7, of the PDR to total emission from 10 -2 to 1. We fix 
the PAH fraction, q pa h, to 4.58, common for solar to sub-solar 
metallicity spiral galaxies (Draine et al. 2007). Errors are computed 
by generating 100 realizations of the photometry using a standard 
Monte Carlo simulation. We find a total dust mass of 5.3 ± 1.0 x 
10 7 Mq, which is factor of ~1.3 greater than a modified black- 
body dust mass with /3=2. The distribution of the dust mass surface 
density is similar with that found by the modified blackbody fit; in 
addition the U m i n map is similar to the temperature map using a 
constant j3. Fig. lAll shows both the dust surface density (left) and 
U m i n distribution using the Draine & Li (2007) models. 

We find that an enhanced 7 can fit the photometry in the nu- 
clear region and the arms, although our variance in these fitted pa- 
rameters is larger than the fitted parameters themselves and thus 
also gives a result consistent with requiring no PDR region. U m in 
is highest in the nuclear region (mean value in nucleus: 11 ± 5; 
maximum value in nucleus: 18 ± 8) and enhanced along the spiral 
arms (mean value in spiral arms: 5 ±1). We also find a similar off- 
set between the peaks in Umin and the dust mass as we have found 
between the dust mass and temperature for the modified blackbody 
model fits. 



APPENDIX B: COMPARISON OF SFR TRACERS 

In this work we used an SFR tracer that was derived from an Ha 
map using a 24 fim map for compact source corrections (see §2.3). 
This tracer is developed by implicitly assuming that the diffuse 
24 /im emission is not associated with recent star formation and 
thus should be removed (i.e. we only use the 24 /im emission as 
a compact source correction). We now examine how these results 
compare to using an SFR tracer that is a linear combination of two 
maps and thus includes this diffuse emission. There are many ways 
to trace sites of star formation and here we explore two alternative 
descriptions using a combination of the far-ultraviolet (UV) from 
GALEX (Gil de Paz et al. 2007) and 24 /jm from Engelbracht et al. 
(2005) as described by Leroy et al. (2008) and a linear combination 



/Hc + (0.031)/ 2 4 Mm . (B2) 

We then convert this into an SFR using Equation 2 (see §2.3) 
and convert this into an SFR surface density. 

Fig. lBll provides a direct comparison of the tracers. On the x- 
axis we plot the SFR used throughout the text, namely Ha, cor with 
a compact source correction using 24 fim emission. For simplicity 
we refer to this simply as Ha. On the y-axis we plot the SFR based 
on the linear combination of two tracers with the UV + 24 fim on 
the bottom panel and Ha + 24 fim on the top panel. In both cases 
the linear combination overestimates the SFR at low values of the 
SFR, with the greatest overestimates being for the UV + 24 fim. 

As clearly stated in Kennicutt et al. (2009), while a linear com- 
bination of Ha + 24 fim may be an appropriate description for the 
SFR as a global quantity, for spatially resolved measures, errors 
will be introduced. Due to diffuse dust emission caused by dust 
heating of older stars, some regions will obtain "spuriously high 
SFRs where little or no star formation is taking place" (Kennicutt 
et al. 2009). For the UV + 24 fim tracer, there is also the diffuse UV 
component and the longer timescale for the UV emission (~100 
Myr) in comparison to Ha (~40 Myr) (see Liu et al. 2011 for a 
discussion and comparison of these tracers). 

In Fig. lB2l we plot dust mass surface density and dust tempera- 
ture versus SFR tracer. We also calculate the correlation coefficients 
between the different pairs. In the case of dust surface density we 
find the correlation coefficients to be 0.75, 0.83 and 0.84 between 
the UV + 24 fim, Ha + 24 ^m and the Ha, cor with a compact 
source correction respectively. In other words, the dust surface den- 
sity is most strongly correlated with the SFR used in this paper. In 
the case of the dust temperature we find an opposite ordering of the 
correlation coefficients with 0.74, 0.69 and 0.62 for UV + 24 fim, 
Ha + 24 fim and the Ha, cor with a compact source correction, re- 
spectively. This suggests that the tracers composed of linear com- 
binations of maps are capturing a diffuse heating component since 
they are more strongly correlated with dust temperature rather than 
dust mass surface density, which is coupled to molecular clouds 
and subsequent star formation activity. Thus, the concern that the 
diffuse 24 emission is not tracing recent star formation seems 
to be valid based on these findings. 



1 We adopt the notation of Draine & Li (2007). 
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Figure B2. Dust mass surface density (upper panel) and dust temperature (lower panel) versus the SFR as traced by a linear combination of UV and 24 fim 
(left) and Ha and 24 fim (centre) and Ha jCor with a 24 fim correction for compact sources (right). Each point represents a pixel of 36 arcsec. The correlation 
coefficients are listed in the bottom right of each panel. 
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Figure Bl. Star formation rate as traced by a linear combination of maps 
(UV + 24 fim bottom panel and Ha + 24 fim top panel) versus a star forma- 
tion tracer of Ha. cor with a compact source correction using 24 fim. Each 
point represents a pixel of 36 arcsec. 



